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Abstract 

In the first part of this series of two papers [Y. Jiao, F. H. Stillinger, and S. Torquato, Phys. Rev. 
E 81, 011105 (2010)], we considered the geometrical ambiguity of pair statistics associated with 
point configurations. Here we focus on the analogous problem for heterogeneous media (materi- 
als). Heterogeneous media are ubiquitous in a host of contexts, including composites and granular 
media, biological tissues, ecological patterns and astrophysical structures. The complex structures 
of heterogeneous media are usually characterized via statistical descriptors, such as the n-point 
correlation function Sn- An intricate inverse problem of practical importance is to what extent 
a medium can be reconstructed from the two-point correlation function 5*2 of a target medium. 
Recently, general claims of the uniqueness of reconstructions using 5*2 have been made based on 
numerical studies, which implies that 5*2 suffices to uniquely determine the structure of a medium 
within certain numerical accuracy. In this paper, we provide a systematic approach to characterize 
the geometrical ambiguity of 5*2 for both continuous two-phase heterogeneous media and their digi- 
tized representations in a mathematically precise way. In particular, we derive the exact conditions 
for the case where two distinct media possess identical 5*2, i.e., they form a degenerate pair. The 
degeneracy conditions are given in terms of integral and algebraic equations for continuous media 
and their digitized representations, respectively. By examining these equations and constructing 
their rigorous solutions for specific examples, we conclusively show that in general 52 is indeed 
not sufficient information to uniquely determine the structure of the medium, which is consistent 
with the results of our recent study on heterogeneous media reconstruction [Jiao, Stillinger and 
Torquato, Proc. Nat. Acad. Sci. 106, 17634 (2009)]. The analytical examples include complex 
patterns composed of building blocks bearing the letter "T" and the word "WATER" as well as 
degenerate stacking variants of the densest sphere packing in three dimensions (Barlow films). 
Several numerical examples of degeneracy (e.g., reconstructions of polycrystal microstructures, 
laser-speckle patterns and sphere packings) are also given, which are virtually exact solutions of 
the degeneracy equations. The uniqueness issue of multiphase media reconstructions and addi- 
tional structural information required to characterize heterogeneous media are discussed, including 
two-point quantities that contain topological connectedness information about the phases. 

PACS numbers: 05.20.-y, 61.43.-j 
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I. INTRODUCTION 



Two-phase heterogeneous media (textures) abound in nature and synthetic situations. 
Examples include manufactured heterogeneous materials (e^., composites, porous media 

and colloids) [11-^], geologic media (e.g., rock formations) [l|, |2|, |5[, cellular materials |6[, 

n 

ecological structures (e.g., tree patterns in forests) [2,], cosmological structures (e.g., galaxy- 
distributions) [sl, |9[, and biological media (e.g., animal and plant tissue) 10|. In general, the 
complex microstructures of random media can only be characterized via certain statistical 
descriptors, such as an infinite set of n-point correlation functions Sn {n = 1,2,...) [l|. 
It is well known that the effective physical properties of heterogeneous media, such as the 



conductivity [ll|, elastic moduli [l2 



.luid permeability |l3|, trapping constant 



3 



and 



electromagnetic wave characteristics 15| , can be expressed in terms of weighted functionals 



of Sn- In particular, S'„(xi,X2, . . . ,x„) gives the probability of finding n points positioned 
at xi,X2, . . . ,x„ all in the phase of interest jl|. For statistically homogeneous media which 
are the focus of this paper, Sn is translationally invariant and hence depends only on the 
relative displacements of the positions with respect to some arbitrarily chosen origin, say 



xi, i.e., S'„(xi, X2, . . . , x„) = S'„(xi2, X13, . . . , xi„) with x^ 



X,-. In such cases, the one- 



point correlation function Si is just equal to the volume fraction of the phase of interest. If 
the medium is also statistically isotropic, the two-point correlation function 6*2 is a radial 
function, i.e., 5*2 (X12) = S'2(|xi2|), providing information about the distribution of pair- 
separation distances. 

Although the Sn can be represented analytically for certain models pj and bounded for 
general media 17|, it is usually not possible to compute all of the Sn in the infinite set. 
Thus, it is desirable to understand the extent to which one can characterize the structure 
and properties of heterogeneous media using lower-order correlation functions, such as the 
two-point correlation function 5*2. A powerful means to study this problem and related 
questions is to employ inverse techniques 18| whereby one attempts to reconstruct (or 
construct) realizations of heterogeneous media that match limited structural information of 
those media in the form of lower-order correlation functions, obtained either experimentally 
or from theoretical considerations 19|, l20|]. In particular, the reconstruction of digitized 
representations of heterogeneous materials from a prescribed two-point correlation function 



5*2 has been receiving considerable attention [2l|-[30|]. An effective reconstruction procedure 





(a) 



(b) 



FIG. 1: (color online). A degenerate pair of two-phase continuous media based on the "kite- 



trapezoid" example given in Ref. 



34l |. The longest distance in the "kite" and "trapezoid" is 



symmetrically placed on the large circle diameter, as adapted from Ref. 



3§. Note that this 



example can be interpreted as two-dimensional (disks) or three-dimensional (spheres). 

enables one to generate accurate renditions of the medium at will and subsequent analysis 
can be performed on the reconstruction to obtain desired macroscopic properties of the 
medium non-destructively. 
It has been well established 
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28 



3lM33| that 5*2 is generally devoid of crucial in- 



formation to uniquely determine the structure of the medium. In particular, in Ref. [33 1 
we employed inverse "reconstruction" techniques to probe the information content of the 
widest class of different types of two-point functions. This set of functions includes the stan- 
dard two-point correlation function S y, s urf ace- void F,^ and surface-surface Fgg correlation 



functions 



35| . lineal-path function L J36|], chord- length probability density function p |37l |. 

. By numerically 



the pore-size function F [l|, and the two-point cluster function C2 
reconstructing two-phase heterogeneous media from these correlation functions, we unam- 
biguously showed that S2 does not suffice for a unique reconstruction and that incorporating 
C2, which is sensitive to topological connectedness information, can lead to a much more 
accurate rendition of the target medium. It is worth noting that in most numerical recon- 
structions, 5*2 of the reconstructed medium only matches that of the target medium within 
certain numerical accuracy [39]. However, it appears that the idea that 5*2 is generally de- 
void of crucial structural information still has not been widely appreciated and claims of 
uniqueness of reconstructions using 5*2 have been made based on numerical studies 22|, l30 |. 



Thus, our work has imphcations for the fundamental problem of determining the necessary 



conditions that realizable two-point functions must possess 40 1 



In the first part of this series of two papers [3J] (henceforth referred to as Part I), we 
have shown via concrete examples the existence of distinct point configurations possessing 
identical sets of pair-separation distances. By appropriately decorating the point configura- 
tions, distinct two-phase media with identical ^2 can be constructed (see Fig. [1]). In general, 
there are many ways to decorate a point configuration, which enables one to obtain a wide 
spectrum of distinct heterogeneous media. In this paper, we focus on characterizing the 
structural ambiguity of pair statistics as embodied in the two-point correlation function 5*2 
associated with individual realizations of statistically homogeneous two-phase heterogeneous 
media in a mathematically precise way, which complements previous numerical studies. In 
other words, we provide a mathematical formulation for heterogeneous-media reconstruc- 
tions and rigorously show that distinct media exist that possess exactly identical two-point 
correlation functions. This objective is quite different from the approximate matching of 5*2 
that one encounters in numerical simulations, which are also discussed. 

The reconstruction of random textures also has a wide spectrum of applications in vision 



research 32|, |4l|, |42|. It is important to distinguish a popular two-point statistics used 
in vision research, called the "dipole histogram" , from the two-point correlation function of 
statistically homogeneous media. The dipole histogram is a particular statistics of the vector 
pair displacements in a finite texture |42|. In such a case, the distinguishable vector pair 
displacement with the largest length can always be identified together with the positions of 
the two points that give rise to the displacement. Similarly, one can then identify the second 
largest displacement and the positions of points possessing this displacement, and so on. 
Finally, the texture can be completely reconstructed from the dipole histogram. However, 
we emphasize that this is possible because the information contained in the dipole histogram 
of a finite texture is not averaged out as that contained in the two-point correlation function 
of a statistically homogeneous medium. In particular, due to the requirement of translational 
invariance, periodic boundary conditions are always imposed for realizations of statistically 
homogeneous media. Thus, only the probabilities of the occurrence of displacements are 
given by 5*2, obtained by averaging over all available positions in the medium, instead of the 
indication of the actual occurrence of particular displacements at certain positions provided 
by the dipole histogram. It is such additional information contained in the dipole histogram 



that enables a unique reconstruction. We will see in the following sections that the position 
average (i.e., volume average) necessarily leads to the structural ambiguity of 5*2. 
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FIG. 2: Illustrations of one-dimensional two-phase random media that form trivial degenerate 
pairs that we do not count. Note that the arrows indicate that the transformations (e.g., mirror 
reflection, translation, phase inversion and these combinations) can occur in either the clockwise 
or counter-clockwise directions. 

To provide a quantitative characterization of the structural ambiguity of 5*2, we first give 
the relevant definitions here. Two heterogeneous media are identical if they possess identical 
sets of n-point correlation functions Sn for n = 1,2,.... Two heterogeneous media are 5*^- 
distinct if they process distinct n-point correlation functions for all n> k. A heterogeneous 
medium is /-fold degenerate if there exist additional (/ — I) media such that all the / media are 
mutually Ss-distinct but possess the same two-point correlation function 5*2. For statistically 
homogeneous media, the above definition of structural degeneracy rules out the possibility 
that two degenerate media are trivially connected by a translation, rotation, mirror reflection 
or any of these combinations (see Fig. [2]) . This is consistent with our definition of degeneracy 
for point configuration in the first part of this series two papers [3J]. Moreover, we consider 
that two media do not form a degenerate pair if they possess phase-inversion symmetry at 
phase volume fractions 0.5, i.e., the morphology of one phase is statistically identical to that 
of the other phase when the two phases are interchanged [l| . 

As pointed out in Part I, given pair statistics in the form of 5*2 that are associated with 
degenerate media, it is impossible even in principle to uniquely reconstruct a medium from 



such 5*2. Furthermore, it is important to determine under what conditions the media possess 
degenerate counterparts. In this paper, we provide a systematic approach to characterize 
the structural degeneracies associated with the two-point correlation function for both con- 
tinuous two-phase heterogeneous media and their digitized representations. In particular, 
we derive the exact conditions for the existence of degeneracy in terms of integral equations 
for continuous media and algebraic equations for the digitized representations. We show 
by examining the derived equations that for statistically homogeneous and isotropic media. 



structural degeneracies generally exist. This explains the long o 
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3served non-uniqueness of 



28l 



29|. We also provide a 



reconstructions using radial averaged 5*2 in many studies [19- 
variety of concrete examples of degenerate two-phase media, including both analytical con- 
structions and numerical simulations, which respectively are exact and approximate solutions 
of the degeneracy equations. These examples include analytically constructed patterns com- 
posed of building blocks bearing the letter "T" and the word "WATER" , degenerate Barlow 
films (i.e., stacking variants of the densest sphere packing in three dimensions) as well as 
numerical reconstructions of polycrystal microstructures, laser-speckle patterns and sphere 
packings. Moreover, we discuss the degeneracy of multiphase media and additional struc- 
tural information that can be incorporated to significantly reduce structural degeneracy of 



heterogeneous media 



33|. 



The rest of the paper is organized as follows: In Sec. II, we introduce the mathematical 
model of two-phase heterogeneous media and derive the conditions of degeneracy. In Sec. Ill, 
we provide the examples of degeneracies for both continuous media and their digitized rep- 
resentations. In Sec. IV, we discuss the degeneracy of multiphase media and additional 
structural information that would be used to reduce degeneracy. In Sec. V, we make con- 
cluding remarks. 

II. CONDITIONS OF DEGENERACY 

In this section, we derive the exact mathematical conditions under which a statistically 
homogeneous medium possesses degenerate counterparts. For a continuous medium, these 
conditions take the form of a set of integral equations of the indicator function (defined 
below) of the phase of interest. For the digitized representations of the medium, which 
have been extensively investigated in various numerical reconstruction studies, the integral 
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equations reduce to algebraic equations. 



A. Continuous Media 



Consider a statistically homogeneous medium M occupying the region V in the d- 
dimensional Euclidean space W'- {d = 1,2,3) which is partitioned into two disjoint phases 
Pj: phase 1, a region Vi of volume fraction (pi and phase 2, a region V2 of volume fraction 
02- It's obvious that Vi U V2 = V and Vi fl V2 = 0. The indicator function X*^*)(x) of phase 
i is given by 

x'"w4' ''"^" (1) 

for i=l,2 with Vj U Vj = V and 

xW(x)+X(2)(x) = l. (2) 

The n-point correlation function Sn for phase i is defined as follows: 

5«(xi,X2,...,x„,) = (X«(xi)X«(x2)...X«(x„)>, (3) 

where the angular brackets "(...)" denote ensemble averaging over independent realizations 
of the medium. The two-point correlation function 5*2 for phase i is defined by 



4^)(xi,X2) = (X«(xi)X«(x2)>. (4) 

As pointed out in Sec. I, for a statistically homogeneous medium, 52 is a function of the 
relative displacements of point pairs, i.e.. 



4^Hxi,X2) = 5?)(x2-xi) = 5«(r), (5) 

where r = X2 — xi. In the infinite volume limit, if the medium is also ergodic the ensemble 
average is equivalent to the volume average, i.e., 

SPir) = \im ^ f 1^^ (x) X(') (x + r) rfx (6) 

y^oo V Jy 



If the medium is also statistically isotropic, S'2 is a radial function, depending on the 
separation distances of point pairs only, i.e., 

4^)(xi,X2) = 5«(|r|) = 5«(r). (7) 

Readers are referred to Ref . [l| for a detailed discussion of S'2 and other higher order Sn ■ 
Henceforth, we will drop the superscript i in S'2 for simplicity. Without further elaboration, 
S'2 is always the two-point correlation function of the phase of interest. 

Now consider a change in the geometry of the region of the phase of interest, say phase 
1, that keeps the volume fraction 0i invariant. The medium after the change is denoted by 
M' and recall that the original medium is denoted by M. The indicator function X' of phase 
1 in M' is then given by 

x'(x)J^ ''^^:- (8) 

[0 XGV{, 
where V( is the region where phase 1 is found in M'. We define the change of indicator 
function as: 

5X(x) =X'(x)-X(x), (9) 

where X is the indicator function of phase 1 in M given by Eq. ([T]) with i = 1. It is obviously 
that 6X can only take the value of 0, -1 and 1. And the requirement that M and M' 
form a degenerate pair with the same two-point correlation function as well as the same 
phase volume fractions will lead to constraining equations on 6X, which are the degeneracy 
conditions we are seeking. 

For simplicity, we introduce a short notation for the integral in Eq. ([6]), i.e., 

S2 = ^X®X. (10) 

The requirement that S'2 is invariant under the change 61 leads to 



82 = 82 = -X' ®X' = -X®X, (11) 



Substituting <^ into ([II]) yields 
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^{I + 6I)^{I + 6I) = ^I®I. (12) 

From Eq. ( 1121) we can obtain that 

I + -51 j ®5X + 5X®(x+ -5X\ = 0. (13) 

Note that given the indicator function X of the original medium M, Eq. (!T3|) should be 
satisfied by 61 for all r within the range of interest. In other words, flT3l) specifies a continuous 
spectrum of equations, one for each 6X at a particular r. Another requirement for degeneracy 
is that the phase volume fractions are conserved, which leads to 

iyJ^X'(x)rfx=lyJ^X(x)dx. (14) 

Substituting ([9]) into ( !H|) . we can obtain that 

(5X(x)rfx = 0. (15) 



de = (17) 



'V 

If the medium M is also statistically isotropic, 5*2 depends only on pair-separation distances. 



S2{r) = ^j^S2ir)de (16) 

where the integral over (i-dimensional solid angle G is to average over all directions of r with 
the same length r and Q is the total solid angle. Thus from Eq. flT3|) . we obtain that 

( X + -6x] ®5X + 5X®lx+ -5X 

Note that ( TT7|) should be satisfied by 5X for all r within the range of interest. 

Thus, it is clear that for the statistically homogeneous medium M , its structural degen- 
eracy M' exists only if Eqs. ( [T3|) and ( !T5|) possess nontrivial solutions; and the degeneracy 
exists only if Eqs. flT7|) and flTSj) possess nontrivial solutions if M is also statistically isotropic. 
In principle, from the solutions of Eqs. (TT3l) [or ( TT71) ] and ( TT5l) together with the indicator 
function X of the original medium, one can obtain the indicator function X' of the degenerate 
medium by applying ([9]). Henceforth, we will call Eq. flTSl) the feasibility condition, since it 
requires that 61 must be feasible in the sense that the phase volume fractions are conserved. 
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Similarly, we will call Eqs. (fT3|) and (fT7|) the invariance condition, since they give a set of 
61 that leave 5*2 invariant under the change. 

B. Digitized Media 

For a digitized medium M in W^, the indicator function takes the form of a finite d- 
dimensional array I = [/a;i...xd] of linear size N, where Ixi...xa = or 1 indicating the phase 
of the pixel at (xi . . . Xd) and x, = 1, . . . , A^ for i = 1, . . . ,d. Consider a perturbation of the 
geometry of phase 1 in the digitized medium M that preserves the phase volume fractions 
(e.g., by exchanging pixels of different phases). The perturbed medium with the indicator 
function I' = [/^^ j.^] is denoted by M'. The change of the indicator function is simply 

61 = 1-1= [4,,,,^ - 4,....J = [6h,...x,], (18) 

which is also a o?- dimensional matrix with linear size N. The entries 6Ixj^...xj^ can only take 
the values —1,0 and 1. 

For statistically homogeneous digitized media, periodic boundary conditions are imposed 
by the requirement of translational invariance. Thus, the indices of the indicator functions 
are modulated by A^, i.e., if (xi + ri) > N the index should take the value (xi+Vi — N) instead, 
where Xi is the original index and r, is the translation displacement along Xi direction. On 
the other hand, if (xj + r^) < the index should take the value (xj + r « + A^) instead. The 
integral equations of the feasibility condition (fTSjl reduce to algebraic equations, i.e.. 



E ^^-i-. = 0' (19) 

Xl,...,Xa 

where the sum is running through Xj = 1, . . . ,N for i = 1, . . . ,d. The equations for the 
invariance condition reduce to 



/ . ^^xi...xa + 2^ SI^x-,+ri)...ixa+ra) + 2^ '^4i...x/A^i+ri)...(xd+rd) - 0, (20) 

(xi,...,Xd)eHi (xi,...,xa)eH2 xi,...,Xd 

where Vi {i = 1, . . . ,d) are integers satisfying |rj| < N/2 and the sets of pixel positions Hi 
and H2 are given by 
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Hi = {(Xi, ...,Xd) I /(a:i+ri)...(xd+rd) = l} , 

(21) 

H2 = {{Xl, ...,Xd) I Ixi...xa = !}• 

If the medium is also statistically isotropic, following Eq. fll7p the invariance conditions 
can be obtained by averaging Eq. (120|) over equivalent directions, i.e., 

/ , I 2^ SIxi...xa + 2^ SI(xi+ri)...{xi+ra) + 2_^ SIxi...xa^I{xi+ri)...{xd+rd) I = 0- 

{ri,...,ra)en \{xi,...,xa)eHi {x^,...,xa)&H2 xi,...,xa J 

(22) 
where the sets of displacements with the same length Q and the sets of pixel positions Hi 
and H2 are given by 

Q = {(ri, . . . , Trf) I rj H h r^ = r^, r G Z, r < iV/2} , 

^1 = {(n, • • • ,rd) I /(a;i+n)...(x,+r,) = l} , (23) 

H2 = {(ri,...,rd) I 4,,,,^, = 1}. 

Thus, similar to their continuous counterparts, two statistically homogeneous digitized 
media M and M' form a degenerate pair if the feasibility condition Eq. ( !T9|) and the in- 
variance condition Eq. fl20|) hold for some non-trivial 51. If the media are also statistically 
isotropic, the degeneracy exists if Eqs. (fT9|) and (!22|) possess non-trivial solutions. Once 51 
is obtained, the degenerate medium M' can be constructed by applying Eq. flTSj) . 

Note that in Eq. ( I22l) . because the number of equations (i.e., the total number of integers 
that are smaller than N'^/4) is in general much smaller than the number of unknowns (i.e., 
the number of the entries Ixi...xa which is equal to A^"'), the algebraic equations fl22l) possess 
multiple solutions. For such isotropic media, it is highly probable to obtain degeneracies of 
the original medium in a reconstruction. This explains the long observed non-uniqueness 
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28 



29|. The sit- 



issue in the reconstruction of heterogeneous media using radial 5*2 [19|, 
nation is more subtle for general statistically homogeneous media, for which the degeneracy 
condition Eq. ( I20l) contains the same number of equations and unknowns. Thus, the chances 
for finding a non-trivial solution of 61 is much smaller but one still cannot rule out the pos- 
sibility of degeneracy. It is clear that the solutions also depend on the original medium, i.e., 

13 



the values of Ixi...xa- In the following section, we will show the existence of degeneracies for 
both continuous and digitized media via concrete examples, which are solutions of the above 
degeneracy equations. 

III. EXAMPLES OF DEGENERACY 

In this section, we will provide concrete examples of structural degeneracy. These exam- 
ples include both exact and approximate solutions of the equations for degeneracy derived 
in the previous section. For continuous and simple digitized media, analytical constructions 
are given. For more complicated digitized media, numerical simulations are employed to 
find degeneracies. For a clear illustration, we will mainly focus on two-dimensional exam- 
ples here. However, the general construction methods can be also applied to find degenerate 
media in any space dimension. 

A. Analytical Constructions 

1. Vector- Argument <S'2(r) 



tWt 




(a) (b) 

FIG. 3: A simple degeneracy example composed of replications of letter "T". 

As discussed in Sec. II. A, the two-point correlation function 5*2 (r) contains information 
on the relative displacements (vector distances) between any two points in the phase of inter- 
est. Consider a two-phase medium, in which one of the phases is composed of replications 
of certain substructures. If the substructures do not possess central inversion symmetry, 
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they can be rotated in a way such that all the relative displacements of point pairs within 
and among the substructures remain the same yet the overall structure (the replication of 
the rotated substructures) is different from the original one. A simple example originally 



proposed in Ref. 4l| is shown in Fig. 121 The substructure, i.e., the letter "T", has been 
arranged into a "pyramid" . By rotating each "T" about an axis perpendicular to the plane 
by 180 degrees, a distinct structure can be produced. Since the relative orientation of the 
position of any pair of letters "T" has not changed, the relative displacements of the point 
pairs, one from each of different "T" also remain the same. It is obvious that the relative 
displacements of point pairs within any single "T" are not affected. Thus, the two distinct 
structures possess identical statistics of relative displacements of point pairs, i.e., identical 
two-point correlation functions 5*2 (r), so they form a degenerate pair. 

Fig. m shows a more sophisticated example based on the same construction rule. The 
substructure now is the letters "WATER" arranged in an "S" -shaped auxiliary box (shown 
in blue), which is then replicated in a complicated way. Note that though the auxiliary box 
possesses central inversion symmetry, the substructure composed of the letters does not. By 
rotating each substructure about the inversion symmetry axis of the corresponding auxiliary 
box, a degenerate structure can be constructed. Note these examples are exact solutions of 
the Eqs. ([15]) and ([13]). 



2. Radially Averaged S2{r) 

Radially averaged 5*2 (r) only contains information on separation distances between the 
point pairs in the phase of interest. In the best case, the complete set of pair-separation 
distances can be inferred from 5*2 (r). Such distance sets have been studied in detail in Part 
I of this series of two papers, where we have shown that a variety of classes of degeneracies 
can be identified that are compatible with the given distance set of a point configuration. In 
particular, we pointed out in Part I that by decorating the degenerate point configurations 
properly, degenerate two-phase media can be constructed. One such example is to decorate 
the 30 degenerate tetrahedra discussed in Part I, e.g., one can place the centroids of congruent 
spheres at the vertices of the tetrahedra and choose the "sphere" phase as the phase of 
interest and the space exterior to spheres as "void" phase. Then the 30 heterogeneous 
media associated these tetrahedra are degenerate. 
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(b) 



FIG. 4: A more sophisticated degeneracy example composed of replications of the letters "WA- 
TER". One "S" -shaped auxiliary box is shown in blue (or light gray in the print version). 

Another example constructed by decorating a point configuration is shown in Fig. |5l 
where the two circular-disk packings form a degenerate pair. The packings contain three 
types of particles: A large number of the circular disks (shown in black) form a cluster 
with central symmetry. Three particles (shown in blue) arranged on a straight line form a 
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(a) 



(b) 



FIG. 5: (color online). A degenerate pair of circular-disk packings in two dimensions. 



substructure possessing mirror reflection symmetry, i.e., the two most separated particles are 
located symmetrically about the one in the middle. The remaining two (shown in red) are 
placed along a line parallel to that of the three blue particles and going through the center 
of the symmetric cluster, such that the set of relative distances between any two points in 
the "particle" phase in the two packings are identical while the two resulting structures are 
distinct. The constructed degenerate pair is an exact solution of the Eqs. (TT5|) and ([1] 





(a) 



(b) 



FIG. 6: A degenerate pair of digitized media based on the same construction rules for the degenerate 
circular-disk packings. 

Note that the aforementioned degeneracy construction rules are also applicable in the 
case of digitized media. For example, applying the construction rule for the degenerate 
circular-disk packings in Fig. |5] to a two-dimensional square-lattice of pixels results in the 
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degenerate digitized media shown in Fig. |6l The only difference is an additional requirement 
that for digitized media the values of the distances must be square roots of integers and the 
symmetry groups of the structures are discrete, which reduces the total number of possible 
degeneracies. 

3. Degenerate Barlow Films 

Another class of interesting degeneracy examples associated with radially averaged 5*2 
involves the Barlow packings, i.e., random stackings of infinite triangular-lattice layers of 
spheres. In a Barlow packing, for a particular layer (i.e., layer "A") there are only two 
possible positions for the next layer, which we denote by the usual code as "B" and "C", 
respectively. Thus, a Barlow packing can be represented by a sequence of "A" , "B" and "C" , 
e.g., ". . .ABCBCA. . .". A Barlow film is a Barlow packing consisting of a finite number of 
layers. The distribution of pair distances in such a geometry depends not only on how many 
layers are present, but also upon whether any chosen pair of layers is the same species ( "A A" , 
"BB", or "CC"), or different species ("AB", "AC", or "BC"). This motivates counting 
"same" versus "different" species for each separation distance between layer pairs. The 
distance distribution (i.e., 5*2) is completely determined by the number of same species at 
all separation distances between layer pairs. 

For an n-layer Barlow film, there are nominally 2" possible configurations. However, 
not all of these 2" configurations are distinct. A basic property of Barlow films is that the 
distance distribution is invariant under a "renaming" operation. For example, the sequences 
"ABCABC" , "BCABCA" and "CABCAB" possess the same structure. In light of this 
"renaming" invariance, we can require that the first two layers for all Barlow films to be 
"AB. . ." without loss of generality. Thus, the number of distinct configurations associated 
with an ra-layer Barlow film that possess the same vertical distances between identical layers 
(i.e., the same distance distributions) is reduced to 2^""^-'. As mentioned in Sec. I, reflection 
(mirror image) symmetry should also be excluded for degenerate pairs. 

We then set out to identify degenerate Barlow film configurations up to n = 15 layers 
by explicitly generating all Barlow films and compare their layer separation distance distri- 
butions. For n < 6, there are no degenerate configurations. An interesting case occurs at 
n = 6, with an inequivalent pair of Barlow films that have the same distance distribution. 
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TABLE I: The number of the same species ( "AA" , "BB" or "CC" ) at different layer-separation 
distances for the 6-layer degenerate Barlow pair "ABABAC" and "ABABCB". (Il is the layer 
separation distance. 
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TABLE IL The number of degenerate Barlow film configurations for n £ [6, 15]. n is the number 
of layers in the Barlow film. Npair is the number of degenerate pairs and Nt^ipiet is the number of 
degenerate triplets. 

Number of Layers Npair NtHpiet 
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and thus the same 5*2 associated with the "sphere" phase. They can be represented by 
"ABABAC" and "ABABCB". As shown in Table [H the two configurations possess the same 
number of same layer species at all layer-separation distances. In addition, they are not 
related by "renaming" and/or reflection operations. 

The numbers of degenerate conflgurations for n G [6, 15] are given in Table Ull We 
note that the flrst degenerate triplet occurs at n = 10, which includes the conflgurations 
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"ABABCBACAC", "ABABCBCBAC" and "ABACBCBCAC". We don't find any four-fold 
degenerate Barlow films for the n values we examined. It can be seen from the table that 
except for the trivial cases with n = 6 and 7 the numbers of degeneracies possess the following 
trends: i) for both the even-layer and the odd-layer films, the number and complexity of 
degeneracies increase monotonically as the number of layers increases; (ii) the number of 
degeneracies associated with an even layer is larger than that of the next odd layer. The 
latter trend is probably due to the odd layers having relatively fewer degrees of freedom to 
exclude the "renaming" and refiection symmetry. 

It is also worth noting that if a honeycomb crystal packing of spheres are used for the lay- 
ers instead of triangular-lattice packings, one can construct degenerate low-density jammed 



sphere packings [43|]. In particular, a honeycomb crystal layer can be obtained by removing 
one third of the spheres from the triangular-lattice layer. These removed spheres are also 
arranged on a triangular lattice, but with lattice vectors that are twice in magnitude of those 
of the original triangular lattice. To obtain honeycomb crystal films from two degenerate 
Barlow films, exactly the same pair distances are removed associated with the removing of 
the spheres. Thus, the two constructed honeycomb crystal films also form a degenerate pair. 
So far, we have not found any degeneracy in the periodic Barlow films. 

B. Numerical Examples 

In many practical applications, one encounters the issue of reconstructing the structures of 
materials from radially-averaged two-point function 5*2 (r), which might be the only available 
information at hand and it is always associated with an inevitable uncertainty. Therefore, 
it is also important to explore the structural degeneracy of 5*2 (r) subject to small uncertain- 
ties. In this section, we employ the Yeong-Torquato heterogeneous material reconstruction 
procedure [19|, 120| to obtain numerical-degeneracy examples of statistically homogeneous 
and isotropic media, which are approximate solutions of the invariance equations fl22l) de- 
rived in Sec. II. B. In particular, the "energy" or mean squared error E for the 



reconstruction problem which has been defined previously 19|, 
given by 



20|, 



m 



33| is 



E = Y,[S2{r) - S,{r)]\ (24) 
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where 5*2 and 5*2 are the target and the reconstructed correlation functions, 
respectively. In the Yeong-Torquato procedure, a stochastic optimization tech- 
nique is used to make E in principle zero. However, in practice E is always a 
very small number, which corresponds to the situation where the uncertainty 
associated with 5*2 is very small. 





(a) 



(b) 
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FIG. 7: (color online). The "polycrystal" microstructures. (a) The target medium, in which the 
volume fraction of the black phase is 0.273. (b) The reconstructed medium using the Yeong- 
Torquato procedure 19| . (c) The two-point correlation functions of the target and reconstructed 
media. The mean squared error E defined by Eq. (j24p is on the order of 10~ . 



We first consider an idealized microstructure of a two-dimensional "polycrystal" shown in 
Fig. [7^. The "polycrystal" contains nine square "grains" , each associated with a particular 
orientation as determined by the direction of the black "stripes" . The two-point correlation 
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function of the black phase is shown in Fig. ^. It can be seen that the oscillations of 5*2 are 
the clear manifestation of the strong spatial correlations between the "stripes" . The recon- 
structed structure is shown in Fig. [7b, whose 5*2 closely matches the target one with a very 
small error, i.e., sum of squared difference between the two, on the order of 10^^. Though 
several stripe-like substructures have been roughly reproduced and a few weakly preferred 
orientations of the substructures can be identified, the reconstruction is distinctly different 
from the target structure in that it contains much more compact substructures and the grain 
boundaries are completely missing. The correlations due to the alternating "stripes" in the 
target medium are mimicked by the correlations between the compact substructures with 
inter-spaces that are approximately equal to those between the "stripes" . 

Another interesting example is the reconstruction of a laser-speckle pattern, which con- 
tains various shaped sub-structures on a broad range of length scales |29| (see Fig. [8^). The 
two-point correlation function (shown in Fig. [8t) clearly does not reflect the multi-scale 
nature of speckle pattern, which monotonically decreases to its long range value very fast. 
The reconstructed pattern is shown in Fig. [8b, with a small error between the correlation 
function on the order of 10"^. The substructures on different length scales in the target 
medium (e.g., the disconnected compact clusters, stripe-like structures and individual pix- 
els) are all mixed up in the reconstruction (e.g., the percolated clusters with different sizes 
and shapes) to reproduce the short-ranged correlations conveyed in 5*2. The reconstructed 
pattern mimics that of the Debye random medium, a famous model structure containing 

nn 

"clusters of all sizes and shapes" with a exponential two-point correlation function [JJ, I19| . 

The above two example clearly show that 5*2 (r) is not able to resolve the details of the 
substructures in a complex medium. But even for simpler structures, 5*2 (r) may still be 
insufficient. Consider a realization of a three-dimensional equilibrium hard-sphere packing 



33| (see Fig. [9^), generated via standard Monte Carlo simulations. The volume fraction 
of the sphere phase is close to the percolation threshold. The reconstruction is shown in 
Fig. [9b, in which the "sphere" phase forms a complex percolating structure. The error 
between the target and reconstructed 5'2(r) is on the order of 10^^^. It is clear that 5*2 
grossly overestimates the percolation of the target phases, which is due to its insensitivity 
to the topological connectedness information of the media. This insufficiency of 5*2 has was 



pointed out previously |l9l. [20 1. 



Note that in the above examples, there are small but non-zero errors between the target 
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FIG. 8: (color online). The laser-speckle patterns, (a) The target medium, in which the volume 
fraction of the red phase is 0.639. (b) The reconstructed medium using the Yeong-Torquato proce- 
dure 191]. (c) The two-point correlation functions for the red phase of the target and reconstructed 
media. The mean squared error E defined by Eq. (j24p is on the order of 10~^. 



and reconstructed 5*2 (r). This is an algorithmic implementation issue and can be removed 
if one used integers instead of floating point numbers in the computer program. However, 
in real- world applications, the obtained 5*2 (r) data would in general suffer small but finite 
errors, no matter how carefully the measuring experiments might be carried out. These 
examples clearly reveal the insufficiency of 5*2 in characterizing multi-scale structures and 
those near percolation. Possible candidate correlation functions to overcome these short- 
comings are discussed in Sec. IV. B. The examples given here are also approximate solutions 



23 





(a) 



(b) 



0.5 



0.1 



— target S^ 
. -. recon. S. 




10 20 30 40 

r (pixels) 



(c) 



50 



FIG. 9: (color online). The hard-sphere packings, (a) The target medium, in which the volume 
fraction of the sphere (blue) phase is 0.446. (b) The reconstructed medium using the Yeong- 
Torquato procedure [l9|. (c) The two-point correlation functions of the target and reconstructed 
media. The mean squared error E defined by Eq. (1240 is on the order of 10^^^. 



of the invariance equations (122|) derived in Sec. II. B. 
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IV. DISCUSSION 

A. Generalization to Multiphase Media 

It is worth noting that although the focus of this paper is two-phase heterogeneous media, 
it is straightforward to generahze the discussion to multiphase media. In particular, consider 
a medium composed of p distinct phases, each associated with an indicator function X*^*) (x) 
[i = 1, . . . ,p) which equals unity when x falls in phase i and equals to zero otherwise. There 
are in total p"^ two-point correlation functions 5*2 (i") including p auto-correlations (i.e., 
when i = j and the two points separated by displacement r falling into the same phase) and 
{p"^ — p) cross-correlations (i.e., when i ^ j and the two points falling into different phases 
i and j). However, there are only p{p — l)/2 independent correlation functions since the p 
indicator functions satisfy the equation 

J(i) (x) + X(2) (x) + ■ ■ ■ + X(P) (x) = 1 , (25) 

and thus, there are only (p — 1) independent X'^*^ 

For every one of the p{p— l)/2 independent two-point correlation functions, similar feasi- 
bility and invariance conditions can be derived, which take the form of integral or algebraic 
equations of the variations of the indicator functions for the corresponding phases, as those 
derived in Sec. II. We will not provide the details of such derivations here. In addition, 
the same arguments concerning the non-uniqueness of the solutions apply in the case of 
multiphase media, and thus one cannot rule out the possibility of structural degeneracy. 

Two simple examples of degenerate three-phase media are shown in Figs. [TH] and [TTl 
which are constructed based on the circular-disk packings and their digitized analog given 
in Sec. IV. A, respectively. A subset of the particles (pixels) have been assigned to the third 
phase, while all the pair-separation distances in the media remain the same. Thus, any pair 
of independent two-point correlation functions of the two degenerate media are identical. 

B. Additional Structural Information 

By examining the degeneracy conditions and constructing concrete examples, we have 
established that two-point correlation function is in general not sufficient to uniquely deter- 
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(a) 



(b) 



FIG. 10: A degenerate pair of continuous three-phase media (black, white and gray) constructed 
based on the degenerate circular-disk packings shown in Fig. [5j 





(a) 



(b) 



FIG. 11: A degenerate pair of digitized three-phase media (black, white and gray) constructed 
based on the degenerate media show in Fig. [6j 

mine the structure of a heterogeneous medium and degenerate structures do exist, especially 
for statistically homogeneous and isotropic media. A natural question is what additional 
structural information could be used to reduce the degeneracy. 

It is notoriously difficult to find the complete answer to the above question. However, 
we have shown in Ref. ^] that it would be a fruitful approach to seek the most sensitive 
statistical descriptors among the various two-point correlation functions (e.g., the surface 



functions and cluster- type functions [ll]) instead of using standard higher-order n-point cor- 
relation functions, as such 5*3. The two-point quantities are superior to 5*3 in that their 
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FIG. 12: A circular-disk system containing small clusters, (a) Target medium, the volume frac- 
tion of the "disk" (blue) phase is 0.532. (b) Reconstruction using 5*2 alone, (c) Reconstruction 
incorporating both S2 and C2. 



determination requires less effort and they are sensitive to nontrivial structural information 
that is not explicitly contained in S3, such as the topologically connectedness information. 

In particular, a superior descriptor proposed in Ref. 33| is the two-point cluster function 
C2(r), which gives the probability of finding two points separated by r falling in the same 
cluster of the phase of interest. A reconstruction example of a two-dimensional circular- 
disk system containing small clusters from both 5*2 and C2 are shown in Fig. [121 It can be 
seen clearly that S'2-alone reconstruction grossly overestimates the percolation of the "disk" 
phase, while the reconstruction incorporating C2 successfully reproduces the disks with small 
clusters. 

We note that proper clustering information can also be used to distinguish certain de- 
generacy examples provided in this paper. Consider the degenerate media bearing a "kite" 
and "trapezoid" structure, respectively, shown in Fig. [H which are obtained by decorating 
corresponding degenerate point configurations. If we increase the size of the circular disks 
while keeping their centers fixed, there exists a critical disk diameter at which the three disks 
in the medium bearing the "kite" structure form a connected "triangle" , which is separated 
from the remaining disk; while in the medium bearing the "trapezoid" structure all the disks 
belong to a single linear cluster. Because the connected "triangle" of disks only exists in 
the medium bearing the "kite" structure, the degenerate pair now can be distinguished, i.e., 
the degeneracy is reduced by using the connectedness (i.e., clustering) information. It is 
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clear that in this case the clustering information reflects nontrivial triangle information in 
the system, which is contained in S3. 

V. CONCLUDING REMARKS 

In this paper, we discussed various aspects of the structural degeneracy associated with 
the two-point correlation function S2 of heterogeneous media. Complementary to previous 
studies, here we provide precise mathematical formulations for structural degeneracy and 
rigorously show that distinct media with exactly identical S2 do exist. In particular, we 
derived the exact conditions for the existence of degeneracy in terms of integral equations 
for continuous media and algebraic equations for their digitized representations. By ex- 
amining the equations and constructing their solutions for speciflc examples, we have well 
established that in general 5*2 is not sufficient to uniquely determine the structures of the 
media, contrary to previous claims of the uniqueness of S'2-reconstructions based on numer- 
ical studies. We have also provided a variety of concrete examples of degenerate two-phase 
media, including both analytical constructions and numerical simulations, which are respec- 
tively exact and approximate solutions of the degeneracy equations. These examples include 
analytically constructed patterns composed of building blocks bearing the letter "T" and 
the word "WATER" , degenerate Barlow films as well as numerical reconstructions of poly- 
crystal microstructures, laser-speckle patterns and sphere packings. It is clearly seen from 
these examples that S2 alone is unable to resolve the details of the microstructures and 
usually overestimates the percolation in the media, which is consistent with the results of 
our recent study on the reconstruction of heterogeneous media using a wide spectrum of 



statistical microstructure descriptors including S2 |33|. The conclusions also apply in the 
case of multiphase media. 

We have pointed out that it is necessary to include additional information to better 
characterize the structure of complex heterogeneous media beyond that contained in 6*2. 
In Sec. IV. B and Ref. [33[, we have shown that the two-point cluster function C2 (i.e., 
topologically connectedness information) is a superior microstructure descriptor for media 
containing compact clusters and its incorporation can significant reduce the degeneracy. For 
more complex media such as those contain multi-scale substructures, structural degeneracy 
associated with higher order S'„ would exist, for example, we may find media that possess 



28 



distinct S4, S5, . . . but identical 5*2 and 5*3. In future work, we will focus on identifying 
such higher-order degeneracies and seeking efficient statistical descriptors (e.g., higher-order 
versions of C2) that can capture the salient features of such media. 
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